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Abstract. We present the spectral and temporal analysis of tlie X-ray flash XRF 01 1030 observed with BeppoSAX. This event 
is characterized by a very long X-ray bursting activity that lasts about 1500 s, one of the longest ever observed by BeppoSAX. 
In particular, a precursor and a late flare are present in the light curve. 

We connect the late X-ray flare observed at about 1300 s to the afterglow emission observed by Chandra and associate it with 
the onset of the afterglow emission in the framework of external shock by a long duration engine activity. We find that the late 
X-ray flare and the broadband afterglow data, including optical and radio measurements, are consistent either with a fireball 
expanding in a wind environment or with a jetted fireball in an ISM. 

Key words, radiation mechanism: non-thermal - gamma rays: burst 



1. Introduction 

Prompt and afterglow emission in GRB show different tempo- 
ral and spectral properties and are usually attributed to different 
mechanisms, more specifically, internal and external shocks, 
respectively. During the prompt emission, the spectrum is hard 
and shows strong spectral evolution from hard to soft. However, 
the afterglow emission is much softer, and its spectral index re- 
mains roughly cons tant during the whole emission from early 
to late observations ( Frontera et al. '2000';' Piro et al. '2002"). The 
transition from one regime to the other takes place from a few 
tens to a few thousand seconds after the burst. In this phase, a 
variety of temporal behaviour is observed in different bursts, 
likely due to the contribution of both prompt and afterglow 
emission. Most intriguing is the presence of X-ray flares, such 
as that observed in XRF 01 1030 . 

Indeed, several bursts observed by BeppoSAX showed the 
presence of X-ray flares fr om tens of second s (e.g., GRB 
970228 iFronteraet all 1998') and GRB 980613 ('Soffitt a et alJ 
[2002,) ) to seve ral minutes (e.g. GRB 01 1 121 and X-Ray Rich 
XRR 011211 llPiroetalJl2005l) ') after the burst. These flares 
have a soft spectrum consistent with that of the late afterglow. 
Furthermore, they connect with the late afterglow emission 
with a power law F oc (f - fo) For early (< 1 minute) X-ray 
flares, the origin of the time to is consistent with the onset of 
the prompt emission. For later (> 10 s) X-ray flares , to needs 
to be shifted to the onset of the flare jPiro et alJ2005l) . 

More recently. X-ray flares taking place in a si milar time 
period have been observed by Swift dGehrels et al...20 04) in 
a larger number of bursts, about one half of the sample 



jO'Brien etalJl2006l) Swift observations are providing a big 
advance in the understanding of this phenomenon and its pos- 
sible relationship with the ce ntral engine. In pa rticular, the dis- 
covery in GRB 050502B ( Bu rrows et al J2005ah of a giant flare, 
~ 700 s after the trigger, with an energy comparable to that of 
the burst itself, suggests that the central engine is undergoing 
long periods of strong activity. 

Swift observations confirm that X-ray flares have a spec- 
trum that is globally softer than the prompt emi ssion, i.e., the 
peak energy Epeak is of the order of a few keV jpalcone et alJ 
'20065. In some cases (GRB 050126, GRB 050219a, and 
GRB 050904), no significant evidence of spectral evolution 
is detected, and the spectrum of the flare is consistent with 
that o bserved in the late afterglow (i Tagliaferri et alJ l2005t 
iGoad et a l. 2006: Burro \ys et al. 2005^). The hght curve con- 
necting the flare with the late afterg low can be reasonabl y 
well-fitted by the (f - fo) " power law ("Tagliaferri et al.'2005'). 
However, in oth er cases (XRF 050406 ( Romano et al . 2006), 
GRB 050502B dBurrows et alJ l2005at Ipalcone et alJ l2006l) . 
GRB 050421, GRB 050607, GRB 050730, and GRB 050724 



llBurrows et al. 2005b)). hard-to-soft spectral evolution was ob- 
served during the flare. 

Several sc enarios were proposed to explain the X-ray flare 
phenomenon JZhang et al. 2006) . Burrows et al. (2005a) pro- 
pose that the central engine releases energy for a long time, 
and internal shocks then produce a long duration prompt emis- 
sion. In the framework of the forward-reverse shock scenario. 
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' We notice that flares and/or re-brightenings are also taking place 
on longer time scales. In the present paper, we focus on flares from 
times of minutes up to ~ 1000 s, i.e., on a time scale similar to that 
observed in XRF 01 1030. 
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iFan & Weil ( l2005l) have shown that, adopting different values 
for the forward and the reverse shock parameters, the reverse 
shock synchrotron radiation can dominate in the X-ray band 
producing a flare. In another scenario, late X-ray flares mark 
the beginning of the afterglow emission, and they are produced 
by a thick shell fire ball (long duratio n engine activity ) through 
an ext ernal shock jPiro et alJl2005h . Very recently, IWu et alJ 
ll2005h have shown that X-ray flares can be produced in the 
context of both late internal and late external shocks. They as- 
sume that the central engine releases energy in two episodes 
(i.e., an early and a late shell are ejected). They applied their 
model to four Swift GRB and found that XRF 050406 and GRB 
050607 flares can be explained both with late internal or exter- 
nal shocks. 

In this paper, we present a complete analysis of XRF 
011030 observed with BeppoSAX. It shows outburst activi- 
ties, with the last detected flare occurring about 1300 s after the 
burst. We investigate the origin of this late X-ray flare. As men- 
tioned above, several models could explain this phenomenon. 
Here we have carried out a detailed analysis of the model in 
which the flare is produced by the interaction of the fireball 
with the external medium. We check if the model can consis- 
tently account for the broadband data - from radio to X-rays. 
We then derive the main parameters of the fireball, including 
the density profile of the surrounding medium. We have also 
tested this model for the X-ray flare occurring in GRB 011121. 

We describe the observations of XRF 01 1030 in Sect.|2and 
perform its temporal and spectral analysis in Sect. l2.1l and l2.2l 
respectively. In Sect. |3] we discuss late flares in the context of 
different variants of the external shock model. In Sect. |3 we 
apply a long duration engine activity (thick shell) model to the 
late flare appearing in XRF 01 1030 and in GRB 01 1 121 , and 
we explain it as the onset of the afterglow emission. In par- 
ticular, in Sect.|4|T] we study the XRF 011030 late afterglow 
emission taking into account the presence of a break occurring 
between 10"* and 10^ s after the burst. In Sect. |5] we use in- 
formation on XRF 01 1030 from the optical and the radio band 
to further constrain the model. Our results and conclusions are 
summarized in Sect.|6l 

2. Observations and data reduction 

The X-ray flash XRF 01 1030 was detected by the BeppoSAX 
Wide Field Camera (WFC) no. 1 on October 30fli, 2001, at 
06:28:02, without any counterpart in the Gamma Ray-Burst 
Monitor (GRBM) ( Gandolfi 2001 ). 

The peak flux is 7.5 x lO^^erg s~'cm~^, and the total fluence 
of the source between 2 and 26 keV is equal to 1.2 x 10"^ erg 
cm"^, consistent with t he typical value o bserved in the same 
range for normal GRB jAmati et al.l2002h . 

The X-ray afterglow of XRF 011030 was identified by 
Chandra in a 47 ks exposure beginning on November 2001, at 
9.73 UT and in a se cond one of 20 ks per formed on November 
2001, at 29.44 UT faarrison et al.l l200lV The locaHzation of 
the X-ra y afterglow was con sistent with the position of a radio 
transient (lTavloretalJ200lh . The radio source was detected on 
November 2001, at 8.80 UT near the centre of the WFC er- 
ror circle at (epoch 2000) R.A.=20:43:32.3, Dec.=-H77: 17: 18.9 



with an error of +1". In this paper, we have used the resul ts of 
the analysis of the Chandra data performed by D' Alessio et alJ 
(I2005ID . The spectrum between 2 and 10 keV is fitted by a power 
law with a photon index F = 1 .72'^^}^ (Table[0. Several optical 
observations were carried out, but none of them succeeded in 
associating an optical counterpart with XRF 01 103 0. The tigh- 
est upper limits are R>21 at 0.3 days after the burst (*Viiav et all 
j^Ol) and R>23.6 at 2.7 days after the burst (Rhoads e talJ 

The precise Chandra localization allowed the association 
of the burst with a faint irregular blue galaxy observed by 
the Hubble Space Telescope and the Keck. The photometric 
observations of this galaxy suggest a redshift smaller than 
z ~3.5, while the low brig htness of the galaxy suggests that 
z > 0.6 teloom et alJ2003^ Since the observations allow us to 
establish only a wide range of redshift values, in our analysis 
we assume z=l. 



2. 1 . Temporal analysis of XRF 011030 

We produced background subtracted light curves of XRF 
01 1030 normalized to the detector effective area exposed to the 
source. The source remained in the field of view (f.o.v.) of the 
WFC for about 1 day. A significant source flux above the back- 
ground level is detected until ~ 1600 s (see Fig. [T). A main 
pulse, lasting ~ 400 s, starts ~ 300 s after a fainter preceding 
event. The main event is also followed by an X-ray flaring ac- 
tivity, 200 s long, which appears ~ 1300 s after the first pulse. 
This X-ray flare has duration and flux similar to the pulse pre- 
ceding the main event. 

As the event is an X-ray flash, its spectrum is soft by def- 
inition. Thus, we cannot expect to find such substantial spec- 
tral differences as those find in GRBs between the phases of 
precursor, prompt emission, and late X-ray flare, i.e., with a 
prec ursor and late fla re markedly softer than the prompt emis- 
sion ( Piro et al. 2005). In any case, due to the similarities of the 
bursting activities that preceded and followed the main pulse 
in XR F 01 1030 to those observed in GRB 011121 jPiro etalJ 
|2005^, in the following we refer to these two pulses as precur- 
sor and flare, respectively. 

After about 1600 s, no signal is detected and we can only 
estimate upper limits on the flux. The light curve of XRF 
01 1030 with the upper limits between 1600 s and 10^ s together 
with the late afterglow emission detected by Chandra is shown 
in Fig. 12 

2.2. Spectral analysis 

We extracted the spectrum between 2-26 keV from the WFC 
data. In the spectral fitting, we tested a simple power law model 
(with and without photoelectric absorption), a broken power 
law model, and a black-body model. The results of our spectral 
analysis are summarized in Tabled All errors are quoted at 1 
cr (68% confidence level). 

The whole spectrum of XRF 01 1030, integrated from s to 
1550 s can be described by a simple power law (Fig.|3j with 
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Fig.l. Light curve of XRF 011030 in the BeppoSAX WFC (2-26 KeV) with a temporal resolution of 50 s. The main pulse 
(300-700 s) is preceded by a fainter event (0-300 s) and is followed by a late X-ray flare (1300-1550 s). 
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Fig. 2. XRF 01 1030 light curve with the BeppoSAX upper lim- 
its on the flux fixed at 3cr above the background fluctuations. 
The late afterglow emission detected by Chandra is shown in 
red squares. [See the electronic edition for a colour version of 
this figure.] 



Q jg, consistent with th e simple power 



a photon index r= 1 . 84+° j ^ 

law photon index F = 1.9 + 0.1 determined bv lHeisel (l2003i 
The fit with this model gives xl - 0.83; 25 degrees of freedom 
(d.o.f.). A fit of the spectrum with a broken power law having 
the photon index Fi fre e to vary and the ph oton index F2 fixed 
to the typical value 2.5 (lAmati et al.ll2 002') did not bring a sig- 
nificant improvement ofx^- Also, the fit with a power law with 
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Fig. 3. vFv spectrum of the total event. The solid line represents 
the fit of XRF 01 1030 with a simple power law model. 



a photoelectric absorption did not bring a significant improve- 
ment ofx^, and led to an upper limit on the absorption column 
density Nh < 1.5 x lO^^cm"^ at z-l. Finally, the fit with a 
black-body model gives a xl value greater than 2, and it can 
thus be rejected. 

We studied the spectral evolution of XRF 01 1030 by divid- 
ing the data into four intervals: the precursor (from 35 s to 280 
s), the first segment of the prompt emission (from 280 s to 500 
s), the second segment of the prompt emission (from 500 s to 
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1200 s), and the flare (from 1300 s to 1550 s). We observed 
only a marginal spectral variability. 

The precursor spectrum is well-described by a simple 
power law with a photon index r=2.61;^Q ^j, marginally steeper 
than the spectrum of the main event. This fit gives ;ifj=0.81; 
25 d.o.f. The precursor can be also described by a black-body 
model with a temperature ^r=0.90;^p [^ keV (see Fig. 13 and 
xl=0.96; 25 d.o.f. This result is interesting because there is 
only one burst, observed b y GINGA, whose sp ectrum is con- 
sistent with a black body (.Murakami et alJl99lh . However, we 
cannot discriminate between these two models. 



3.1. The Fireball Model: the standard "thin" shell case 

In a "standard" approach, the Fireball Model assumes a thin 
shell (Sari & Piran 1999) that expands with spherical symme- 
try either in a constant density medium or in a wind profile en- 
vironment. In this framework, the emitted flux reaches its max- 
imum at the deceleration radius ro and then starts to decrease. 
However, it does this too slowly to be consistent with the onset 
and decay of the flare. This appears clearly in Fig.|5j where the 
calculated light curve for a thin shell fireball expanding in an 
ISM is shown. 




Energy (keV) 



Fig. 4. vFy spectrum of the precursor. The solid line is the fit of 
the precursor with a black-body model. 

The first and the second parts of the prompt emission are 
both fitted by a simple power law; for the first part, the photon 
index is r=1.78^Q and;\fj=1.24; 25 d.o.f., and for the second 
one, r=1.63+°;^|^ and ;rv=l-48; 25 d.o.f. 

Finally, the late X-ray flare is also fitted by a simple power 
law. Its spectrum is marginally steeper than those of the main 
event, as we also found in the case of the precursor: r=2.10;|;y 
with;i'2^1.51; 25 d.o.f. 

3. The late X-ray flare in the context of external 
shock models 

Among the difl'e rent models proposed for X-ray flares 
(IZhanget al."200^), we choose to analyse in detail some of the 
models based upon an external shock origin, motivated by the 
spectral similarity observed in the flare and afterglow phases, 
straightforwardly accounted for in this scenario. A detailed 
analysis is carried out to check the capability of the model to 
account for the whole set of broadband data. 

In what follows, we first try to explain the late flare of XRF 
01 1030 as being due to ext ernal shock in a "stan dard" fireball 
model (i.e., thin shell case. ISari & Piranl (Il999l) '). with a con- 
tinuous or discontinuous density profile. Since the flare cannot 
be described by this model, considering the similarity of XRF 
01 1030 wifli GRB 01 1 121, we finafly explain it by shifting the 
origin of time to to the instant of the flare, which corresponds 
to a thick shell fireball. 



Fig. 5. X-ray light curve of XRF 011030 for a standard thin 
shell fireball expanding in an ISM with £53 = 0.03, Fq - 45, 
n - \, Be - 0.3, eb - 0.05, and p - 2.1. 



Moreover, we note that a thin shell fireball also fails to ex- 
plain the emission at about 1000 s. In such a case, the deceler- 
ation time of the fireball t^ec is greater than the duration of the 
central engine activity te„g, thus prompt and afterglow emission 
are well separated, and one would expect no emission between 
these two phases. 

3.2. Model with a discontinuous density profile 

A discontinuous density profile can be produced by a variable 
activity of wind emission or by interaction of the wind bubble 
with the external uniform medium. When the fireball expands 
in a medium characterized by a sudden increase of density, 
one could expect that a larger number of photons is produced 
and the flux increases quickly. This is true only when the elec- 
tron cooling frequency Vc is greater then the observational fre- 
quency Voij. On the contrary, when Vc < Vots, the emitted flux is 
independent of the density pr ofile both during the slow co oling 
and the fast cooling regimes JPanaitescu & Kumar II2OO0I) . Let 
us then verify if the cooling frequency could be located above 
the X-ray range at the time of the flare. 

A fireball expanding in an IS M decelerates at a time tjec 
dPanaitescu & Kumar I2OOOI) : 



(1) 
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Table 1. Results of the spectral analysis of XRF 011030. The models used are: Power Law (PL), BroKeN power law (BKN), 
Power Law plus photoelectric ABSorption (PL+ABS) and Black Body (BB). 



name 


interval 


model 


photon 


Nh 




kT 


FluX2-26 keV 


xi 


V 




(s) 




index F 




(keV) 


(keV) 


[erg ■ cmT^ ■ .s"'] 






total 


0-1550 


PL 


1.84+"" 

-0.16 












0.83 


25 


event 




BKN 


1.77t'9 

-U.iJ 




< 11 




1.3 ■ 


. 10-9 


0.81 


24 






PL+ABS 


1.88+°-?I 

-0.14 


< 1.5 - lO^"* 






1.2 ■ 


■ 10"' 


0.86 


24 


precursor 


35-280 


PL 


2.6 r"-^? 

-0.61 








J. o 




0.81 


25 






PL-hABS 


2.44+2-06 

-0.44 


< 7 - 10^- 






J.J 




0.87 


24 






BB 








0.90+"!' 

-0.15 


3.7 - 


10-'" 


0.96 


25 


burst 


280-500 


PL 


1 78+"-''' 

^ - '°-0.16 








2.5 ■ 


10-'' 


1.24 


25 


part 1 




BKN 


1 tQ+0.23 

— U.Z/ 




9.1+^-? 




2.3 ■ 


10-' 


1.08 


24 







PL+ABS 


2 23+"--'6 


< 4.4 - 10^^ 







2.1 ■ 


10-' 


1.16 


24 


burst 


500-1200 


PL 


1 63+"-^'^' 








8.9- 


10-'" 


1.48 


25 


part 2 




BKN 


39+1--*'' 




< 4.4 




7.8- 


10-'" 


1.6 


24 






PL+ABS 


1 Qi+0.74 


< 8.6 - 10^- 






8.1 - 


10-'" 


1.65 


24 


X-ray late 


1300-1550 


PL 


2 lO+O-**' 








4.9- 


10-'" 


1.51 


25 


flare 




BKN 


4.18 ± 10 




< 3.0 




5.8- 


10-'" 


1.7 


24 






PL+ABS 


2 10+^'''' 


< 1.3 - 10^3 






4.9- 


10-'" 


1.72 


24 


afterglow' 


(9.24-9.71)- 10^ 


PL+ABS 


1.72 + 0.20 


2 96+" "' - 10^' 






5.8- 


10-'^ 


0.76 


9 



'From lD'Alessio et aljj2005l) : the flux is in the 2-10 keV range. 



During the deceleration phase (for t > t^pr), the cooUn g 
frequency Vc is given by jPanaitescu & Kumar l2000l) : 

Vc ^ 3.4 X lO^^E^l'^-^e'g^ilq^'^ Hz. (2) 

Except for very low values of the Lorentz factor, Fq < 30, 
the flare occurs during the deceleration phase and the cooling 
frequency v^. is given by Eq.|2l This equation indicates that, for 
typical XRF energies and parameter values, at the time of the 
flare, f3 - 1, the cooling frequency Vf could be higher than the 
observational frequency Vobs (that is, in the X-ray band) only 
for small values of density, n < 0.07. During the deceleration 
phase, Vc decreases with time and will pass below the X-ray 
range at later times. 

If the fireball starts its expansion in a wind density profile, the 
deceleration time tdec is tPanaitescu & Kumar -2000.) : 

tdec - 6.67£5iA;,L2ro,l s (3) 
and during the deceleration regime 
Vc ^ 3.77 X 10"'£^(^A;22e^^/2ji^2 

Also in this case, for A, < 10"^, v^ can be higher than Vobs 
at the flare time. Moreover, now Vc increases with time, and 
therefore the X-ray afterglow emission will remain sensitive to 
density variations. Thus, for a suitable range of parameters, at 
the time of the flare the X-ray emission can be sensitive to den- 
sity. However, the duration and amplitude of the flare are not 
co nsistent with the k inematic upper limit recently established 
bv lloka e t al.' (2005) on the flares produced by the interaction 
of the fireball with density discontinuities. In particular, if we 
assume to observe GRB on axis, this upper limit is ( loka et alJ 
l2005l) : 

^<V_L^, (5) 

Fy 5 vFv t - to 



where ' ~ {vilvc)'''' ^^^^ is the fraction of cooUng energy 
and F/vFy ~ (vobs/yc)^''^^^^^ for v^ < Vobs- The flare duration 
is At ^ 200 s in XRF 011030, and the time of the flare oc- 
currence is {t - to) ^ 1300 s, where the time is counted from 
fo, i.e., in the case of a thin shell, from the initial trigger. Thus 
At/t ~0.15, and with typical parameters values, Eq.|5limplies 
AFy/Fy < 0.25, while from the X-ray data of XRF 01 1030 we 
derive AFy/Fy ~ 3.6. 

We point out that the case discussed above applies 
only to a density disco ntinuity with a shell geometry. 
iDermer & Mitmanl ( 1 19991) have shown that a clumpy medium 
would be able to produce high variable light curves through ex- 
ternal shock if the clouds radius is very small in comparison to 
their distance from the central engine. This process can explain 
X-ray flares up to thousand of seconds ( iDermer 120051) . 

3.3. Long duration engine activity: the thicl< shell 
model 

In the following, we show how we can describe the flare in the 

context of the external shock scenario by shifting the origin of 

the time to to the onset of the flare. From a theoretical point of 

view, the onset of the external shock depends on the dynamical 

regime of t he fireball that is str ictly related to the "thickness" 
I 7" — I 1 

of the shell (Sari & Piran 1999), i.e., to the duration of the en- 
gine activity. In fact, a shell is defined as being thin or thick 
depending on its thickness A - cte„g, where t^ng is the dura- 
tion of the engine activity, and also on its initial Lorentz factor 
Fq. The shell is defi ned to be thick if (£/nm,,c2)i/3r^8/3 < ^ 
dSari & Piran 11199"^ . For our purpose, we rewrite the above 
equation s ubstituting the deceleration time given in Eq. (3) of 
IPiro et all t2005) . obtaining t^ng > tdec for the thick shell con- 
dition. Most of the energy is transfeiTed to the surrounding 
material at tdec for thin shells or at teng for thick shells. In the 
latter case, the peak of the afterglow emission therefore coin- 
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cides with te„g. Also, the afterglow decay will be described by 
a power-law only if the time is measured starting from the time 
at which the inner engine turns off, i.e., fo - tf„g. According to 
iLazzati & Beeelman .1.2006) . this should happen when the cen- 
tral engine releases most of the energy during the last phase of 
its activity. In this context, the flare would thus be produced by 
the external shock caused by an energy injection lasting until 
the time of the flare occurrence, i.e., teng ~ 1300s. The hypoth- 
esis of external shock for the flare offers a straightforward ex- 
planation of the spectral similarity with the late afterglow data. 
We also notice that in this model the early after glow emission 
is mix ed with the (internal shock) GRB emission dSari & Piran I 
Il999l) . In this context, the emission observed at 1000 s can be 
attributed to internal shock while the flare represents the onset 
of the afterglow emission. 

To develop our model, we used the prescriptions of the 
so-called standard fireball model (Panaitescu & Kumar 2000). 
They offer analytic solutions only at distances greater than 
the deceleration radius ro. But we are also interested in the 
so-ca lled coasting and transition phases (Zhang & Meszaros 
l2004l) because we want to study and to reproduce the shape of 
the flare, its rise included. As mentioned above, we have taken 
into account the thick shell variant by introducing a time shift 
to, i.e., implying that most of the energy for the external shock 
is carried at fo ~ tg„g. We thus numerically solved the basic 
equations of the fireball model. The program requires the pa- 
rameters of the model, namely the initial value of the Lorentz 
factor of the relativistic shell Fq, the energy value in unity of 
10^^ erg £53, the electron population index p, the fraction of 
energy going into relativistic electrons Se, the fraction of en- 
ergy going in magnetic field sb, and the density of the external 
medium n (cm^^*) or A» (cm"') as input. The density profile is 
described by the law n = 3.0 ■ lO-'^A.r"', where in the case of 
an ISM, s -0, while in the case of a wind profile environment. 

When not stated otherwise, we have taken £53=0.03, as- 
suming that all the kinetic energy is converted in y-rays and 
that the redshift is z = 1 . From our spectral analysis, we find 
p - 2.1. To determine the other parameters values, we per- 
formed a study devoted to understanding how they influence 
the calculated light curve. 

We investigated the effect of model parameters on the X- 
ray light curve produced by a thick shell fireball with spherical 
symmetry expanding in an ISM. The origin of the time, fo, is 
shifted to the instant of the flare, 1300 s. We show the X-ray 
flux between 2-10 keV obtained by numerical integration of the 
specific energy flux. 

First we find, according to the fact that the X-ray emis- 
sion is typically above the cooling frequency, particularly at 
late times, that the density n and sb only have a marginal effect 
on the normalization of the X-ray light curve. Figure|6lshows, 
for example, the effects of the density n of the external medium 
in which the fireball expands. Differences are appreciable only 
at early times, when the observational frequency Vobs is smaller 
than the cooling frequency Vc- 

The effects of the parameters £53 and are presented in 
Fig. The normalization of the X-ray light curve depends 




Siooo 10'' 10'' 10'' 10^ 



Time (s) 

Fig. 6. Effects of density n of the external medium on the X- 
ray light curve with the origin of the time shifted to fo=1300 s. 
For the red curve n=5 cm^^, for the blue curve n-l cmT^, and 
for the green curve «=0.1 cmT^ . The X-ray light curve peak 
shows a small increase with n. The other model parameters are 
£53=0.03, ro=85, e,=0.01, £^=0.05, and p^2.2. 

mostly on the product of £53 and e^., and follows a roughly 
linear dependence (for p x 2.2). 




Time (s) 



Fig. 7. Effects of the energy £53 and Sg on the X-ray light curve 
with the origin of the time shifted to fo=1300 s. For the green 
curve £53 = 0.03 and Se - O-l, for the blue curve £53 =0.3 and 
Ee = 0.01, for the Hght blue curve £53 = 0.03 and Se = 0.01, 
for the red curve £53 = 0.003 and = 0.01, and finally for 
the orange curve £53 = 0.03 and Se - 0.001. The other model 
parameters are n=l, ro=85, 6^=0.05, and p-2.2. Note, as the 
emitted flux increases with £53 and Sg, that these parameters 
have about the same amount of influence in determining the 
normalization factor of the X-ray light curve. 

Figure |8l shows the effects of the initial Lorentz factor To. 
Its value influences both the height and the wideness of the 
peak in the X-ray light curve. The greater To is, the higher and 
narrower the peak is. For t » fo, i.e., when the deceleration 
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phase has been reached for the different values of Fq, the light 
curve is independent of Fq. 




^000 10* 10- lo'' 



Time (s) 

Fig. 8. Effects of the initial Lorentz factor Fq on the X-ray light 
curve with the origin of the time shifted to fo=1300 s. The blue 
curve was obtained for Fq - 200, the green curve for Fq = 85, 
and the red curve for Fq - 30. Note, as the peak in the curve 
increases with Fq, the other model parameters are £53=0.03, 
„=1, e,=0.01, eB=0.05, and p^2.2. 



4. External shock from long duration engine 
activity in XRF 011030 and GRB 011121 

We applied the model described in the previous section to XRF 
Oil 030 and also to GRB 011121. In the case of XRF 01 1030, 
we studied the event both in an ISM and in a wind profile en- 
vironment, producing the calculated light curves and finding 
a family of solutions corresponding to several choices of the 
model parameters. In Figs.|9land^|we report two possible so- 
lutions for a fireball expanding in an ISM and in a wind profile 
environment, respectively, with the origin of the time shifted to 
the onset of the flare, to = 1300 s. SmaU changes (< 10%) of to 
do not appreciably modify the results. 

We find that the calculated light curves can describe the 
flare, both for a fireball expanding in a wind profile environ- 
ment and for a fireball interacting with a uniform medium. 
These two light curves do not fit the late afterglow data; this 
will be discussed in the next section. 

In the case of the flare of GRB 011121, we computed the 
light curve only f or a fireb all interacting with a wind den- 
sity profile. In fact. lPiro etal. ( 2005) estabhshed that the GRB 
011121 X-ray and optical data are consistent with a fireball 
expanding in a wind environment due to the temporal decay 
observed in these two bands. Using their parameters and shift- 
ing to to the onset of the flare, we find the light curve of Fig. II 1! 
The model describes the flare and the late afterglow, in agree- 
ment with the analysis made bv .Piro et al.. (.2005.) . 




100 :o()() 10* lo"' 10^ lo'' 

Time (sec) 



Fig. 9. X-ray light curve of XRF 011030 in an ISM obtained 
shifting the origin of the time to the onset of the flare, fo - 
1300s. The model parameters are £53 = 0.03, Fo = 1 10, n = 1, 
Be = 0.2, SB - 0.05, and p -2.1. 




100 1000 lo* 10° 10'' 10'' 

Time (s) 



Fig. 10. X-ray light curve of XRF 01 1030 in a wind obtained 
shifting the origin of the time to the onset of the flare, to - 
1300.S. The model parameters are £53 = 0.03, Fo - 50, A, - 
0.05, = 0.1, = 0.05, and p = 2.1. 



4.1. The interpretation of tlie breai< in the light curve of 
XRF 011030 

In Fig. 12 we note that the backward extrapolation of the late 
afterglow flux detected by Chandra is not compatible with the 
upper limits observed by BeppoSAX, suggesting the presence 
of a temporal break. 

First we considered the possibility that the temporal break 
is related to a spectral break, i.e., to the passage of the cooling 
frequency Vc in the X-ray band. We first studied the case of a 
wind density profile, when Vc increases with the time as f'^^. 
Initially, Vc can be smaller than v„bs, but there will be an instant 
at which it becomes greater than Vobs- This marks a break in the 
light curve, which becomes steeper by 6a = 0.25. 

The observational data of XRF 011030 suggest that the 
break occurs between lO'* and 10^ s after the burst. We thus 
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Fig. 11. X-ray light curve of GRB 011121 for a fireball ex- 
panding in a wind with the origin of the time shifted to the 
instant of the flare, fo - 250 s. The model parameters are 
£53 = 0.28, To = 130, A. = 0.003, Se = 0.01, eb = 0.5, 
and p - 2.5. 

require that Vc passes in the X-ray band during this temporal 
range. Eq. (0} for a given time of the break Tt links A» with 
Eb- We derive Eg and To values to reproduce the light curve of 
the flare as described in the previous section. A» is constrained 
from late X-ray data and also optical and radio data (see next 
section). Finally, the corresponding value of eb is constrained 
from Eq. We first attempt to find a broadband solution with 
the isotropic energy fixed to £53 = 0.03 (see Sect. 13. 3> . In this 
case, we have problems fitting the radio data. We find a set 
of model parameters able to describe the emission observed in 
the X-ray and optical bands, but the corresponding radio light 
curve is always below the observational data. When the fireball 
expands in a stellar wind, the flux in the radio band goes as: 

FyOcEl'-^^A, E^J^_y (6) 

The parameters Ee and eb are determined by X-ray and 
optical data (see the next Sect, for more detail). Then, if 
we keep £53 = 0.03 to obtain the right normalization of 
the radio light curve, we need to increase the wind den- 
sity. On the other hand, this will also cause the normaliza- 
tion of the X-ray and optical light curves to increase and 
surpass the observational data. This has motivated us to as- 
sume an efliciency 77 - 0.1 to convert the kinetic energy re- 
leased by the central engi ne in y-rays. This choice i s sup- 
ported b y several authors. Guetta, Spada & Waxman ' ("2001") 
and Kobavashi & Sari (200 1 ) have argued that internal shocks 
convert the energy with an efficiency 77 ~ 0.1 - 0.5, and this 
was also recently supported by Swift observations (Zhang et al. 
I2OO6; Granot. Konigl & Piran 2006 ). Under this assumption 
(i.e., £53 = 0.3) we find that it is possible to explain the ra- 
dio data jointly with the X-ray and optical data, obtaining a 
broadband solution. In Fig. [21 we show the calculated X-ray 
light curve. 

It is interesting to note that before the break (vc < Vobs), i.e., 
during the flare, the expected photon index isr = (;?/2H-l) = 




100 1000 10* 10^ 10^ lo'' 

Time (s) 



Fig. 12. X-ray light curve of XRF 01 1030 in a wind with the 
origin of the time shifted to fo = 1300s. This solution also takes 
into account the radio and optical data. The model parameters 
are £53 = 0.3, Fo = 60, A, = 0.055, e, = 0.02, sb = 0.001, and 
p = 2.1. We have assumed the efficiency in the conversion of 
the kinetic energy to be 77=0. 1 . 

2.05, consistent with the value found in our spectral analysis 
(Sect.Q. After the break (v^ > Vobs), the expected photon index 
is F = [(p - l)/2 + 1] = 1.55, which agrees with the value 
F = 1 .72 + 0.20 found bv lP'Alessio et all ilOO^ in the analysis 
of the Chandra data. 

With regard to the ffreball temporal evolution after the break, 
Vf. > Vghs ( Panaitescu & Kumar -20021) : 

F oc t-^'P-'^'\ (7) 

and we expect that the temporal decay index is a2 - 1.325. 
The analysis of Ch andra data sho ws that after the break a2 - 
2.25 ±0.60 (D'Ales sio et al.l2005l) . these values are marginally 
consistent. 

Similar considerations can be made for a ffreball expanding 
in an ISM. In this case, the cooling frequency Vc decreases with 
the time as r^^^. Supposing Vc > Vgbs before the ffare, there will 
be an instant when becomes smaller than Vobs and a break 
occurs. 

Now, after the break, the temp oral decay is slower than the 
case of a wind density proffle dPanaitescu & Kumar l2002i) : 

F oc r(^''-2)/4 

and the temporal decay index expected after the break is a2 - 
1 .075, not consistent with the Chandra data. Moreover the spec- 
trum after the break steepens, in disagreement with the spectral 
data. 

In the ISM case, it is therefore even more difficult than in 
the wind case to explain the late afterglow emission without 
introducing a jet structure. The emission coming from a rela- 
tivistic shell with jet symmetry is similar to the one of a spher- 
ical ffreball, as long as the observer is on the jet axis, and the 
jet Lorentz factor 7 is greater than the inverse of its angular 
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spread 00 (lRhoadslll997l) . During its expansion, the fireball 
collects a growing amount of matter; thus, the Lorentz factor 
■y decreases and there is an instant when y < 0^'. At this time, 
the sideways spread of the jet becomes important and the ob- 
served area grows more quickly. This leads the flux to decrease 
more rapidly whit respe ct to the s pherical c ase, and we exp ect a 
break in the light curve dSari et al..l999.) . .Sari et alJ ( ll999) cal- 
culated that at high frequencies the flux decreases like f both 
when Vohs > Vc and when Vobs < Vc- Thus, with the electron 
population index p - 2.1, the predicted temporal behaviour 
agrees with the two Chandra observations. Once the sideways 
expansion of the jet become s important, the cooling fr equency 
Vc is constant with time (P anaitescu & Kumar 11200 l l) and the 
spectrum should not evolve. In our data, the XRF 01 1030 spec- 
tral evolution is only marginally significant; in fact, the photon 
index F of the power law fitting the flare is consistent within 
the errors with the photon index of the power law describing 
the afterglow (TableQ). We therefore carried out a comparison 
between the model and broadband data, i.e., taking into account 
the optical and the radio information discussed in the next sec- 
tion. 

In this case, we find a solution that nicely describes all the 
data without requiring an efficiency r] in the conversion of the 
kinetic energy in y-rays (see Fig.^jfor the X-ray light curve). 
We notice that even if the jet model has an additional free pa- 
rameter with respect to the spherical fireball model for a jetted 
fireball, the model parameters are, still well-constrained. This 
is mostly due to the passage of the cooling frequency Vc in the 
X-ray band. At the start of the observation Vc > Vobs, at about 
10^ s, the cooling frequency becomes smaller than Vobs- After 
this instant the X-ray and optical flux follow the same law, and 
this well constrains the model parameters. The constraints on 
the model parameters given by optical and radio information 
are discussed with more detail in the next section. 




100 1000 10* 10^ lo" 10^ 

Time (s) 



Fig. 13. X-ray light curve for a jetted fireball expanding in an 
ISM with the origin of the time shifted to fo = 1300,?. The 
model parameters are £53 = 0.03, Fq = 130, n - 'b, Se — 0.29, 
SB = 8 • 10"^ p = 2.1, and T/, = 8 • 10^ i. 



5. Broadband analysis of XRF 011030 afterglow 
data 

No optical counterpart has been detected for XRF 011030. 
From among all the optical obse rvati ons, we considered th ose 
performed by'Viiav et al.' feOOl") and'Rhoads et al.' feOOll) be- 
cause they are t he most constraining. The upper limits are 
R>21 llViiavet alJ[2001l) and R>23.6 dRhoads et alJl20ol at 
0.3 and 2.7 days after the burst, respectively. We correct mag- 
nitudes for the reddening due to the absorption of our Galaxy, 
finding R > 20.4 and R > 23.1 (corresponding to an optical 
flux Fy^opn = 1.79 X 10"^** erg cm"^ s"' Hz"' and Fy^„p,2 = 
1 .69 X 10~^^ erg cm~^ s~' Hz"') for the two observations. In the 
radio band, iTavlor et al.l 1I2OOII) associated a transient source 
with a flux = 1.81 x lO"^^ erg cm"^ s"' Hz"' about 10.5 
days after the burst with XRF 01 1030. 

For a jetted fireball expanding in an ISM, a solution that 
accounts for X-ray (Fig.ll3>. optical (Fig.O, and radio (Fig. 
[15} is given by £53 = 0.03, Fq = 130, n = 5, = 0.29, 
SB = 8 ■ 10 ^ p = 2.1, and T/, = 8 ■ 10^ s. We show the op- 
tical and radio light curves corresponding to this set of model 
parameters in Figs, and respectively. We investigated 
how well the parameters are constrained, with particular re- 
gard to the density. The density is mostly constrained by the 
data below the cooling frequency, in this case optical and ra- 
dio, and whether times are greater than about lO'* s (that is the 
time when a spectral break occurs), also X-ray data. After 10"* 
s, the emitted flux in the X- ray and optical band is given by 
dPanaitescu & Kumar I2OO0I) : 

This same relation applies for the radio data because, at 
the time of the observation, the injection frequency v, is below 
or very near the observational frequency Vobs- Thus, the nor- 
malization of the light curve in one of the three observational 
bands also determines the normalization of the light curve in 
the other two bands. This causes the model parameters to be 
well-constrained. 

In the case of the spherical fireball expanding in a wind, 
the break has to be self-consistently described (i.e., without the 
addition of a free parameter). Consequently, the model parame- 
ters are also well-constrained, £53 - 0.3, Fq = 60, A» - 0.055, 
Sg = 0.02, sb = 0.001, and p = 2.1. The corresponding X-ray, 
optical, and radio light curves are shown in Figs. \^ and 

ini 

6. Summary and conclusions 

In this paper, we have presented the temporal and spectral anal- 
ysis of the X-ray flash XRF 01 1030 observed by BeppoSAX. 
This event is o ne of the longest in the BeppoSAX sample 
din't Zand et al.k2002) . with a duration of ~ 1500 s. In particu- 
lar, along with the main pulse, we find a precursor event and a 
late X-ray flare. 

While the spectrum of the main burst is not consistent 
with a black body, we cannot exclude this model for the pre- 
cursor This result could be due to the lower statistics avail- 
able, but it is nonetheless interesting to note that, so far, there 
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Fig. 14. Optical light curve of a jetted fireball expanding in an 
ISM with the origin of the time shifted to fo = 1300 s. The 
model parameters are £53 = 0.03, To = 130, « = 5, = 0.29, 
fiB = 8 • 10"^ p = 2.1, and T/, = 8 ■ 10^ s. 
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Time (s) 

Fig. 15. Radio light curve of a jetted fireball expanding in an 
ISM with the origin of the time shifted to fo = 1300 s. The 
model parameters are £53 = 0.03, To = 130, « = 5, = 0.29, 
SB = 8 ■ 10-^ p = 2.1, and T;, = 8 ■ 10^ s. 

was on ly one example of a p recursor consistent with such a 
model jMurakami et aljfiggil) . This feature could be associ- 
ated with the cocoon formed by the jet emerging at the sur- 
face of the collapsing massi ve star (.Ramirez-Ruiz et al...2002.: 
IWaxman & Meszaros l200:^. 

After the launch of the Swift satellite. X-ray flares appear 
to be a common feature in GRBs light curves. This has favored 
the development of a large number of models to explain X- 
ray flares' origin, both in internal and external shock scenarios 
jZhang et al. 2006 ). X-ray flares observed by BeppoSAX in 
GRB 01 1 121 and in XRR 011211 have spectra similar to that 
of the late afterglow. This similarity can be straightforwardly 
accounted for in the framework of the external shock, i.e., the 
flare represents the onset of the afterglow and it is connected 
with the late afterglow emission with a power law. For GRB 
011121 and XRR 01 121 1, a connection with the late afterglow 
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Time (s) 

Fig. 16. Optical light curve of a spherical fireball expanding 
in a wind with the origin of the time shifted to fo = 1300s. 
The model parameters are £53 = 0.3, To = 60, A, - 0.055, 
Be - 0.02, bb - 0.001, and p - 2.1. We have assumed the 
efficiency in the conversion of the kinetic energy to be rj-Q. 1 . 
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Fig. 17. Radio light curve of a spherical fireball expanding in 
a wind with the origin of the time shifted to fo = 1300s. The 
model parameters are £53 = 0.3, Fo = 60, A, = 0.055, Be = 
0.02, bb - 0.001, and p -2.1. We have assumed the efficiency 
in the conversion of the kinetic energy to be rj-0. 1 . 

is acceptable only by shifting the origin of the time fo to the in- 
stant of the flare ( Piro et a l. 2 005) . This implies a long d uration 
engine activity (thick shell fireball. lSari & Piranl ( ll999i) V 

We find that this scenario fits nicely with the observational 
data, including the X-ray flare and the late broadband radio, op- 
tical, and X-ray afterglow observations of XRF 011030 . The 
latter, performed by Chandra, indicate the presence of a tem- 
poral break occurring between 10'* s and 10^ s after the burst. 
We carried out a detailed modelling of the data, finding good 
agreement with observations for a spherical fireball expanding 
in a wind medium and for a jetted fireball expanding in an ISM. 
In the first case, the temporal break is explained by the passage 
of the cooling frequency in the X-ray band. We cannot exclude 
that the flare observed in XRF 01 1030 is due to internal shocks. 
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However, this would require an engine activity that is tuned to 
track the peak energy of the emission during the prompt and 
flaring phases. 

In the context of the external shock scenario, the time shift 
to ~ 1300 s is due to a long lasting central engine activity that 
remains active until the time of the flare, with t he most of the 
energy re leased at the end of the emission phase dSari & Piran I 
ll999H Lazzati & Beaelman 2006). In this case, the peak of the 
afterglow emission coincides with the flare, and the afterglow 
decay will be described by a self-similar solution by counting 
the time from the instant the inner engine turns off', i.e., when 

We also find that a thin shell fireball cannot describe the 
flare in a continuous density profile, in fact, in this case, the flux 
rises and decays too slowly to describe the shape of the flare. 
Also, with a discontinuous density profile, it is very difficult 
to produce flares as high and narrow as that observed in XRF 
011030, unless one assumes that the fireball is expanding in a 
clumpy medium ( D ermer & Mitrn g n ..1999: .Dermer .2005) . 

Recently, Swift observations have shown the presence of 
late X-ray flares in several other events. Some of these flares 
appear to have a spectral behaviour consistent with that of 
late X-ray flares observed by BeppoSAX, i.e., a soft spec- 
trum substantially well-differentiated from the hard, prompt 
emission typically attribute d to in t ernal s hocks (GRB 050126, 
GRB 050219a IxagHaferri et all l2005l) . and GRB 050904 
llBurrows et alJEo05b) ). GRB 050126 and GRB 050219a ap- 
pear to follow a (f - toy power law, with to ^ 100 s, reason- 
ably well. Instead in other flares, e.g XRF 050406 and GRB 
050502B, the hardness ratio su ggests a sp ectral evolution re- 
sembling the prompt emission ( Burrow s et al..2005bl This be- 
haviour has been thus interpreted by Burr ows et alJ 112005 ah as 
due to internal shocks produced by a long duration central en- 
gine activity. However we notice that in some cases it could be 
possible explain hard-to-soft spectral evolution also in the con- 
text of the external shock scenario. When the fireball expands 
in a wind the cooling frequency Vc increases with the time and 
there will be an instant when it enters in the X-ray band. At this 
instant the spectral index changes from /3 - p/2 to/3 = {p-l)/2 
and the spectrum becomes harder of a factor 0.5. 

Swift observations showed also the presence of multiple 
flares in GRB light curves on relatively small time scales, from 
~ 100 s up to ~ 1000 s. For example GRB 050421 show s 
two successive flares within about 150 s i Godete^alJ 2005j). 
GRB 050607 has two flares within 500 s ( Pagani et al. 200^) 
and GR B 050730 has three flares within 800 s ( Bun'ows et aL 
l2005bl) . In these cases a thick shell fireball can be successful 
to explain only one of the flares appearing in the light curve, 
i.e. only one flare can be identified with the beginning of the 
afterglow emission. The other flares can be attributed to inter- 
nal shocks or to the interaction with a clumpy medium in the 
framework of the external shock scenario. In conclusion we 
note that the present data suggest the existence of two cate- 
gories of late X-ray flares differentiated by their spectral be- 
haviour. It is also interesting to note that both in the framework 
of the internal shocks scenario and in that of external shocks, 
the explanation of late X-ray flares requires a central engine 
that remains active about until the time of the flare. 
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